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Abstract 



In a previous paper we demonstrated that Bethe's equations are not suf- 
ficient to specify the eigenvectors of the XXZ model at roots of unity for 
states where the Hamiltonian has degenerate eigenvalues. We here find the 
equations which will complete the specification of the eigenvectors in these 
degenerate cases and present evidence that the sl2 loop algebra symmetry is 
sufficiently powerful to determine that the highest weight of each irreducible 
representation is given by Bethe's ansatz. 



I. INTRODUCTION 



The XXZ model with periodic boundary conditions, defined by 

H = ~\ X>K + i + ojfoj+i - cos 7 a^ +1 ), (1.1) 

Z 3=1 

where ex* is the i Pauli spin matrix at site j and j — L + 1 = 1, has long been studied by 
means of what has come to be called Bethe's equation 

sinh \{ Vj + ij) \ L _ 1 sinh \( Vj - y t + 2ij) 
sinh — ij) J , = 7 sinh ~(vj — v / — 2ij) 
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where S z = | J2f=i &j is a conserved quantum number. The eigenvalues of are given 
by 

£ = ^-2^'^ (1-3) 

2 ~[ cosh Uj — cos 7 

and the eigenvectors are given in terms of the Vj. 

The authors who originally derived these equations!!! concentrated on the energy and 
eigenvector of the ground state and this computation has been done with great rigor and 
clarity However, the derivations seem to superficially extend to all eigenstates and for 
almost three decades there have been extensive effortJffl made to study these excited 
energy states and for generic values of 7 it has been proven0 -i! that the solutions to ( |1.2|) 
do in fact give all the 2 L eigenvalues of ( |1 . 1|) . 

On the other hand it was shown by BetheB in his original study of the case 7 = 
that for the Heisenberg antiferromagnet the equation ( |1.2| ) with finite Vj does not give all 
eigenvalues. Furthermore it was realized!! as early as 1973 that the equation ( |1.2|) will not 
specify solutions of the form where both the numerator and denominator simultaneously 
vanish and that these special solutions do indeed occur in the "root of unity case" 

7o = ^vr (1.4) 

where r and N are relatively prime and 1 < r < N — 1. Therefore there are indeed values 
of 7 which are not generic in the sense of0B where Bethe's equation ( |1.2j ) with finite Vj is 
known not to give all eigenstates of the system. 

Recently it was demonstrated!^ that when the root of unity condition ( |1.4|) holds that 
the Hamiltonian (|1.1| ) has an invariance under the s/2 loop algebra and that this symmetry 
leads to degenerate multiplets of energy eigenvalues. Moreover we subsequently showed^ 
when 7 — > 70 that there are solutions v k to (|1.2|) which have N members of the form 



v k = a + i2k~f with 1 < k < N (1.5) 

where a is in general complex with —tt/N < Ima < tc/N. We referred to these solutions 
as exact complete N strings and in ref.u we numerically presented many examples. These 
solutions give a contribution of zero to the energy independent of a by virtue of the identity 

jv 1 

^ cosh(a + ik2nr/N) - cosm/N 

and thus are responsible for degeneracies in the energy eigenvalue spectrum. 

These exact complete N strings all have the property that in the equation ( |1.2|) they give 
rise to factors 0/0. Therefore the equation (|1.2j ) will not be able to determine the parameter 
a of the exact complete N string ( |1.5| ). It is for this reason that in ref.0 we said that Bethe's 
equation is incomplete at roots of unity. 

However, it is clear that the missing equations can be obtained in principle by setting 

7 = 7o + e (1.7) 
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in ( |1.2|) and carefully taking the limit e — > 0. It is the purpose of this paper to carry out this 
construction which will complete the specification of the Bethe's roots v^. 

For a fixed value of S z > any solution will have 1 < k < ~ — S z . If the solution 
contains n of the complete exact N strings there will be n Q other roots where 

n = | - S z - nN (1.8) 

We denote these roots as v °. We will call these roots "ordinary" roots because we will show 
that they satisfy the Bethe's equation (where we will let L be even) 

/ sinh + yy ) \ L _ ^ sinh ± (tf - vf + 2i-f ) 

{ sinh \ (v) -ho)) ~ t\ sinh \ (u° - vf - 2i 7o ) ' 1 ' ' 

This equation for t>° does not involve the parameters ay of the exact complete N strings 
and by definition does not have any ambiguous terms of the form 0/0. These ordinary roots 
Vj may be real or may have imaginary parts which are organized into strings (see ref.0 for 
many examples). For the purposes of this paper this information is not needed and is thus 
not indicated in the notation. 

The simplest case with complete exact N strings is the state which does not contain any 
ordinary roots. For a chain of length L these special states have S z = |f — nN where n is 
the number of exact complete iV strings. The n parameters a m are determined from the n 
equations 

N-i i 

£ sinh L -(a m + (2k + l)i 7o ) x 
k=o 1 

I j n N-1 y \ 

Lcoth-(a m + (2A; + l)z7 )-2^ £ coth -{a m - or, + i2j lQ )\ =0 (1.10) 

\ 1 3=0 1 

In general when there are ordinary roots in the state the n parameters a m are determined 
from 

N-1 ^ N+2k+l 

Y, sinh L -(a m + (2k + 1)270) I] Piiphn) x 

k=0 ^ l=2k+4 

■y n N-1 ^ 

L coth -(a m + (2k + 1)270) ~ 2 51 H cotn 9 ~ a i + iZjlo) 

1 1=1 j=0 1 

f>oth X -(a m - + 2A»7o) + coth ^(a m - u° + 2(fc + l)«7o))^ = (1.11) 



where 



n a ^ 

P k (a m ) = II sinh n(om - u? + 2^7o) (1.12) 



and we note the periodicity P^A^a/J = P^ct^). 
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In the special case N = 2 the product Y\iPi(oi m ) is replaced by one and setting 

z m = exp(a m ) Cm = exp(>^) (1.13) 



the equation fll.llj ) reduces to 



/ / - 1 i \ ' lit + l ^ L 2 + ( l ) L ( Zm l ) L _ a v z ™ + z t - 9 V z ™ + ^ - - n n i i \ 

H z m + 1 -) /.. IT _ JVL _ T! (i\LI 7 _A\L 4 2^ 72 3 2 2^ _ 2 _ 72 ~ U " (MI) 

-1- £j -t- yt) \z m b) !=1 z m z ; i=1 2 m C,; 



Z = l 



It is to be noted that this equation is not the equation which is obtained by replacing the 
righthand side of ( |1.2| ) by (— 1) L / 2_ I S,Z I _1 which is what would be obtained if we "formally" 
set 7 = tt/2. We also note that the case N = 2 has been treated in ref.il in quite a different 
context. 

These results will be derived in section 2. We will see that in order to determine the 
parameters aj we need to expand (O) to order e 2 and that the equation ( |1 . 1 1| ) is obtained 



as a consistency condition arising from the vanishing of the secular determinant of a set of 
homogeneous linear equations. 

We will discuss the consequences of our result and the possible relations with the evalu- 
ation parameters introduced in ref.0 in sec. 3. In particular we present compelling evidence 
that for N > 3 the sZ 2 loop algebra is powerful enough to determine not only the structure 
of the degenerate multiplets but the highest weight vectors themselves. On the other hand 
the highest weight vectors are already known from the Bethe's ansatz solution of the XXZ 
spin chain. Therefore we conclude that the Bethe's equation (|1.2p is in fact contained in the 
irreducible representations of U(sl?) at level zero. 



II. DERIVATION OF (1.11) 



We begin our derivation by first considering (|1.2| ) directly at 7 = 70 where Vj is an 
ordinary root v® and split the product over / on the right hand side into the n a contributions 
from ordinary roots and the nN contributions from the N strings to find 



'sinh ~(d? + ry 




y sinh I (v^ - 

*l sinh ±(y° - v? + 2i lQ ) " sinh |(j - a m - 2i{k - l) 7o ) 

l\ sinh \{v? - v? - 2* 7o ) ii x W sinh - a m - 2i{k + l) 7o ) ' 1 ' ' 



The product over k gives unity independent of a m and thus (|2.1| ) reduces to ( |1.9|) as desired. 
In order to derive ( |1.11| ) we cannot directly set 7 = 70 in (|1.2| ) and let vj = Vj ik = 



ctj + 2iA;7o because factors of 0/0 will occur on the right hand side. Therefore we first write 
sinh L ^{v jjk - i"i)q{v j)k + 227) + sinh L ^(v j:k + i^)q{v j;k - 2ry) = (2.2) 

with 
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n j n TV j 

g(uj >fc ± 2ij) = JJ sinh - - v° ± 2ry) JJ JJ sinh -(u j)fc - ^ m ± 2ry). (2.3) 

i=l ^ 2=1 m=l ^ 

We then define for the roots which become exact N strings 

v jtk = a 3 + 2ik lQ + evfl + e 2 v$ . (2.4) 
and for all other ordinary roots 

vi = vf + €v} 1] . (2.5) 

We proceed in several steps. 

Expansion of sinh L \ {vj$ ± ij)- 

We expand sinh L \{vj :k ± 17) to order e as 

sinh L ^(uj,fe ± ij) ~ sinh L i {otj + 2i/c7 + ei^ ± 2(70 + e)) 

~ sinh L \ [otj + (2A; ± 1)270] U + ^(^2 ± i) coth ^ + (2k ± l)i 7o A . (2.6) 

Expansion of q(vj jk ± 2ry). 
To expand q(vj >k ± 227) we write 

n n a 

q{v jtk ± 217) = f(v j>k ± 217) JJ gi(v j<k ± 217) ]J ^i( u i,fc ± 2i 7) (2-7) 



with 



i=i 



TV x 

/(u j>fc ± 2ry) = [J sinh - v jtm ± 2^7) 

m=l Z 
TV ^ 

gi(v jtk ± 2ry) = J| sinh - v hm ± 2ry) with j ^ I 

m=l 1 

}n(v j>k ± 22 7 ) = sinh -(v j>k - v? ± 22 7 ) (2.8) 



and expand /, g t and h t separately. We find 

TV 

sinh 

2 



1 

f(v j7k ± 217) = [] sinh 9 ( 2iA; 7o + + <?vfl - 2im 7o - euj^ - e 2 vf m ± 2i( 7o + e 

m=l Z 
TV 

= ]~[ (sinh[?(A; — m ± 1)70] + 



m=l 



\(-ir N ^ sinhmno) (e(v§ - v$ ±1 ± 2.) + e 2 (^ - v% ±1 ) 

TV-1 1 \ 

+(«$ - t&i ± 2^e 2 £ - «<JL (1+I) ± 2i) cothiZ 7o J , (2.9) 



N x 

gi(v j>k ± 217) = ]J sinh -(wj,fc - v hm ± 2ry) 

m=l ~ 



I] sinh ~( aj -a t + 2i{k -m± l) 7o + e[u$ - ± 2z] 



m=l 

N 



_^r(jv-&±i) -q ginh 1 ^ . _ ^ + 2im l0 ) x 

m=l ^ 



/ N 1 1 \ 

( 1 + e E 2 - w £ ± 2z i coth 2 - a < + 2 ^ - m ± 1)7o) ) (2 - 10) 



and 



hi(v jtk ± 2* 7 ) = sinh - (otj - v? + 2i(k ± l) 7o + e(vfl - ± 2i 
= sinh - (aj — Vi + 2i(k ± 1)70) x 



1 + <vfl - v\ 1] ± 2i)i coth - «,° + 2^(fc ± l) 7o ) ) • (2.11) 



where we have defined 



(i) _ (i) (i) _ (i) 



(2.12) 



Expansion of equation (2.2) to order e 

We now are able to expand Q2.2|) . Because /(i>j,fc ± 2*7) vanishes as 7 — > 70 the leading 
term is of order e. It is convenient to define 

Zj.fc = vfl - vf k+1 + 2% {-x j , k -x = vfl-vfl_ 1 -2i) (2.13) 

with 

iV 

J2xj,k = 2iN (2.14) 
fe=i 

and also let 

<f> k {aj) = sinh L -(or, + (2k - 1)270) (2.15) 

where we note the periodicity 4>k+N(&k) = <f>k(oc k ) Then we find that for each value of j the 
equation (|2.2| ) reduces to leading order in e to a set of N homogeneous equations for the 
variables Xj tk with 1 < k < N 

(j) k (aj)P k +i(aj)xj jk - fa+xfa^Pk-iia^Xj^-.! = (2.16) 

where we recall the definition of P k (ctj) of (|1.12 ). The determinant of the coefficients of this 
set of linear equations vanishes for all values of the still undetermined aj. Thus equations 
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( [2.16| ) are consistent and with the normalization condition ( |2.14| ) we find Xj jk in terms of ctj 
as 

r .„ = 2lN fa+ijaj) (217) 
hk K(aj) P k ( aj )P k+1 ( aj ) 1 • ] 

where 

fc=i ^feWJ-nfe+ilaiJ 

We of course have not yet found the desired equations to determine ctj. To do this we 
need to expand (|2.2|) to one more order in e. 

Expansion of equation (2.2) to order e 2 

We now define 

yi,k = vfl-vfl +1 (2.19) 

and find 

4>k{aj)Pk+i{aj)y jt k - <j>k+i(aj)Pk-i(aij)yj,k-i = ~ R k (2.20) 

where 

i? fc = 4 1} + 4 2) + 4 3) + 4 4) (2.2i) 

with 



R k ] = j(4fc ~ l ) coth \i a i + ( 2k ~ 1 )*7o]0fe(aj)-Pfe+i(a J >i ) fc 



^( v i,fc + *) coth \\ a 3 + ( 2fc + l)ilo)<f>k+i(aj)Pk-i{®j)xj,k-i, (2.22) 



1 Ar ~ 1 

^1 2) = tE ( u S - + 2i ) cothi/7o0 fc (a j )P fe+ i(a j )x iifc 



4 
1 



i=i 

AT-1 



7 51 ( W S ~ ~~ 2i) cothi77o0 fc+ i(a J -)i , fc-i(aj) a: j,fc-i> ( 2 -23) 



R * = 7 Z) XI - W S + 2 «) «>th - «J + 2(fc - m + l)i7 O )0 fc (a i )P fe+ i(a i )x ii fc 

4 (=1 m=l / 

Mi 

i n N ^ 

-tEE ( u i,fc ~ W S - 2 «) coth n(«j - "J + 2(fc - m - l)i7o)0fc + i(a i )P*_i(a j )a:i,fc-i 

4 l=i m =l Z 

(2.24) 
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and 

1 n ° 1 

R k ] = 7 X>jJ - V ^ + 2l ) COth o " V l + 2 ( fc + WloWkMPk+lMHk 

4 J=l z 

1 n ° 1 

- I X>j,fe - v i - 2l ) coth o("i - v ° + 2 ( fc - lJ^^fc+iKO^A-iKO^fc-i- (2-25) 
4 z=i ^ 

The left hand side of ( |2.20| ) is identical with the left hand side of the order e equation 
(with Xj^ —* Hj,k) and thus the N x N determinant of the coefficients of yj^ vanishes. This 
implies that the Rk must satisfy the following constraint 

E M<*j h 1 P & j \ \ = 0. (2.26) 
The consistency equation for atj 

It remains to substitute the expression for Rk ( |2.21|) and Xj^ (|2.17|) into (|2.26|) . We 
consider the four terms R^ separately. 
For R^ we find 

fc= i 4 t^i V 2 fc+ i 

-KJ + i) coth -a, + (2fc + 1>7„ J — — (2.27) 

2 <Pk J 

NL 2i /, m s , l r , , v , 
= ~K(^ S - coth 2^ + ( 2fc ~ ^ 



+ i) coth -[a,- + (2A; + l)i T o]J (2.28) 
iiVX ^ 1 

= "oF? — V E coth -[a, + (2k + l)i T o]%* (2-29) 

Zii \ a 3> k=l L 

= E coth - a, + (2fc + l)i 7o — , 2.30 

where to obtain ( |2.29| ) we let k — > A; + 1 in the first term of ( |2.28| ) and then use the definition 
of Xj t k and to obtain ( |2.30| ) we use (|2.17| ). We remark that even though vfy appears in 
R^(aj) of (|2.22|) only the differences Xj^ appear in (|2.30| ). 



For R^? we find 



- («$ - - 20 cothzho ^^^ ^ fc -i) (2.31) 

■jy N N-l 

77777-- r E E " wjJ+i-H + 4 c ° th «7o (2-32) 
fc=i z=i 

(2.33) 
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where to obtain ( |2.32|) we have used ( |2.17|) and to obtain ( |2.33|) we have let k — > k + 2(1 + 1) 
in v ( i 1 l_ l _ 1 and have used the antisymmetry of coth 2/70 under I —+ N — I. 



For Rjp we find 



» R^(a 3 )P k (a 3 ) 
k= i 0fc(«i)0fc+i(«i) 



+ 2z) coth — a i + 2i(k — m + 1)70) 



1 AT n N 

tEEE 

^ fc=l (=1 m=l \ 

■ - t{m - 2i) coth \((Xj — ai + 2i(k — m— 1)70) 



Pk{aj)Pk+i(aj 

<t>k+l{0ij) 



Pk-i{aj)Pk(aj 
4>k(°<-j) 



~ X j,k-l 



(2.34) 



iN 



2K(a 



N n N 

EEE 

3) k=l 1=1 m=l 



y ffc ~~ "fi + 2 coth - (a,- - ai + 2i(k - m + 1)70) 



( v fl ~ v \]i - 2i ) coth ~ a i + 2i ( k ~ m ~ l)7o)) 



2iV 



2 n AT 



1 



w ^EE coth 9 ( a J - «J + 2im 7o ). 
A 1=1 m.=i z 



(2.35) 
(2.36) 



We note that even though vf\ appears in R^(ctj) it has canceled in the final expression 



Finally for P^ we have 



^ R^(g j )P k (g j ) 
k =i 4>k{aj)4>k+Mj) 



N n a 



E E ( («j3 - «i W + 2 coth - «? + 2 ( fc + l)<7b) 



fc=i i=i 



Pfc(Qj)Pfc + i(qj 
0*+i(a*) 



1 



- ( V # - ^ - 2z) coth -(a 3 - «? + 2(k - l^lo f^f^ x^ 

2 fc iati) 



(2.37) 



•Aj AT n a , 1 

" EE mi - + 2 *) coth o(«i - + 2 ( fc + !)^) 



2 #K, fc =u=i 
- ( V $ - - 2i) coth - «,° + 2{k - l)i 7o )) 

Z '7V - n ° / 1 
-— — £ £ coth - v? + 2(k + l)i 7o ) 

+ coth ^(a, - vf + 2/ci7 )^ 

iN N n , y 

J2J2[™^-(a 3 -v? + 2(k + l)i lo ) 



(2.38) 



(2.39) 



2K(a 



3) k=i 1=1 



+ coth - (a 3 — vf + 2ki'-f ) S J 



2iN 4>k+i(aj) 
K(a 3 ) P fe (o; i )P fc+1 (a i ; 



(2.40) 



where we note that (for which we have not derived an expression) has canceled out in 
going from (gjg ) to (|2~39| ). 
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Thus we may use (ggjj) , (g33[) , ( ^36|) and (|CT|) in ( ^26|) to obtain 



1 * 
V- 

^K) £^i p k(a j )P k+ i(a j/ 



{Lcoth^( aj + (2A; + l)z 7o ) 



n ° / 1 1 
5] coth - (Qj - u,° + 2(& + l)i 7o ) + coth -(a,- - u,° + 2^7 ) ) } 

n N ^ 

2 cotn 9 ( a i ~ a « + 2 ^ m 7o) = 0. 

i=i m =l ^ 



(2.41) 



The result (|1.11|) of the introduction follows immediately. 



III. DISCUSSION 

The equation ( |1.9| ) for the ordinary roots vf and the equation ( |1 . 1 1| ) for a*, replace the 
Bethe's equation ( |i~2"D at roots of unity where undefined factors of 0/0 occurred. Thus for 
the case of 7 = rn/N we have succeeded for the first time in completely specifying the 
parameters Vj which occur in the Bethe's Ansatz wave function as given by Yang and Yang! 
and in the auxiliary matrix Q(v ) in Baxter '^1,0 functional equation for the transfer matrix 
of the six vertex model. We have also numerically verified for the cases N = 2 — 5, r = l 
that the equations ill .9|) and ( |1 . 1 1| ) reproduce the numerical results previously obtained by 



other means in ref.l 

There are several other features of our computation which deserve to be discussed in 
detail. 

First of all even though we have expanded ( |2.2|) to order e 2 and have found an explicit 



solution (|2.17 ) for the variable Xj^ ( 2.13 ) which depends on the differences 



V': 



(1) _ .,(1) 
j,k V j,k + 1 



have not obtained expressions for the first order corrections v^l themselves which will in 
general contain a constant independent of k which is not present in Xjj.. This constant plays 
a role in the first order correction identical to the role which atj plays in the zeroth order 
solution and is determined from the consistency equation needed for the determination of 
the e 3 corrections. This pattern of the consistency equation for order e n+2 being needed to 
completely specify the roots to order e n is a general feature for all orders in e. 

Secondly we point out that the parameters a introduced in ( |1.5| ) as part of the specifica- 
tion of the complete exact iV string do not have to be real and that if a is complex then we 
see from taking the complex conjugate of ( |2.41| ) that a* will also be a a solution (because 



all ordinary roots appear in complex conjugate pairs). These complex values of a are a new 
feature of the solution of Bethe's equation which were first seen in our0 previous numerical 
treatment of the problem where many examples were exhibited. 

Thirdly we acknowledge that to completely solve the eigenvalue degeneracy problem we 
need to be able to classify and count all solutions of ( |1.9| ) and ( |1.11| ). We studied this 



degeneracy in ref. iandi in terms of the finite dimensional representations of the s/2 loop 
algebra and we remarked that from the theory of affine Lie algebrasil it is known that these 
representations are specified by what are called "evaluation" parameters and that these 
parameters are roots of the Drinfeld polynomial. However these evaluation parameters are 
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not the same as the roots ay studied in this paper and we note in particular that for a 
multiplet containing 2 m degenerate eigenvalues there are only m evaluation parameters even 
though there are 2 m solutions for aj all of which have the same n ordinary roots. Thus it 
would seem that there is a sense in which there is more information in the aj than what is 
needed for the solution of the degeneracy problem for the eigenvalues of the Hamiltonian 
(and the transfer matrix). 

To explore further the nature of the information contained in the evaluation parameters 
we have explicitly computed them for a chain of L = 12 and N = 3 and have found that 
the evaluation parameters of each multiplet are different. This means that there are no 
isomorphisms between multiplets and thus we reach the striking conclusion for this case 
(and we believe for all cases with N > 3) that the sfa loop algebra is not only a symmetry 
algebra which classifies states into degenerate multiplets but that it is sufficiently powerful 
to determine all the highest weight vectors as well. Therefore the complete decomposition of 
the representation of the loop s/2 algebra into finite dimensional irreducible representations 
must produce as highest weight vectors the same vectors which are determined by the Bethe's 
equation ( |1.2| ) when used in the Bethe form of the wave function!. This relation between 
finite dimensional irreducible representations of the loop algebra of SI2 and Bethe's ansatz 
has not been previously noticed. 

Finally we mention that unlike the previous discussion^ of the degeneracies of the XXZ 
Hamiltonian at roots of unity in terms of the SI2 loop algebra for which the generalization 
to the XYZ model is unknown, the considerations of this present paper have a very natural 
extension to the XYZ model. Details of this extension will be published elsewhere. 
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